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TECHNICAL NOTE 2786 


EQUIVALENT PLATE THEORY FOR A STRAIGHT MULTICELL WING 
By Stanley U. Benscoter and Richard H. MacNeal 


SIMMARY 


A structural theory is developed for the analysis of thin multi- 
cell wings with straight spars and perpendicular ribs. The analysis 
is intended to be suitable for supersonic wings of low aspect ratio. 
Deflections due to shearing strains are taken into account. The theory 
is expressed entirely in terms of first-order difference equations in 
order that analogous electrical circuits can be readily designed and 
solutions obtained on the Cal -Tech analog computer. 


INTRODUCTION 


In the process of designing thin supersonic wings for minimum 
weight it is found that a convenient construction with aluminum alloy 
consists of a rather thick skin with closely spaced spars and no stringers. 
Such a wing deflects in the manner of a plate rather than as a beam. 
Internal stress distributions may be considerably different from those 
given by beam theory, particularly for torsional or eccentric loading. 

A thin-walled multicell wing may be regarded as an elastic shell 
in which the webs, or wall segments, act as membranes. By introducing 
appropriate simplifying assumptions it is possible to reduce the equa- 
tions for such a shell to the form of plate equations provided that one 
simple criterion is satisfied. The chordwise cross sections must have 
a horizontal- axis of symmetry. This condition is usually satisfied in 
supersonic wings which do not contain stringers. 

The theory which has been developed is of the simplest form that 
could be reasonably expected to give results of an accuracy that would 
be satisfactory for engineering purposes. The loading is assumed to 
consist of a set of vertical concentrated forces acting at the inter- 
sections of the ribs and spars. The distributed load on the wing must 
be replaced, in some rational manner, by an equivalent set of concen- 
trated loads. 

In order that a solution may be obtained on an analog computer, 
the structural theory is expressed entirely in terms of first-order 
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difference equations. The design of the analogous electrical circuit 
is explained and illustrated. 

This paper was prepared at the California Institute of Technology * 

under the sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. 

SYMBOLS 

Ar area of rib ( shear web only) 

A s area of spar .(shear web only) 

a,b,c,d components of displacement which midpoints of a panel take 

when panel is put in shear _ 

C shearing stiffness of a spar or rib 

D bending or twisting stiffness of. an element 

Dij' = XipEl/(l - n 2 ) 

% 

Dji' = XjMEl/(l - H 2 ) * 

E Young ' s modulus 

G shearing modulus of elasticity 

h thickness of wing 

I moment of inertia of skin per unit of width 

Ir moment of inertia of vertical web of rib 

I s moment of inertia of vertical web of spar 

i number of a rib 

number, of a spar 

increase in a function across ith rib and jth spar, 
respectively 


J 

Ai,Aj 
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bending moment 

number of a bay along a spar; also a mass 
number of a bay along a rib 

increase in a function across mth bay of a spar and nth bay 
of a rib, respectively 

concentrated load 

shear flow 

flexibility 

thickness of skin 

twisting moment 

shear 

deflection 
spanwise coordinate 
chordwise coordinate 
vertical coordinate 

angles center lines of a panel rotate through when panel Is 
put in shear 

rotation of a normal 

shearing strain 

length of a bay 

Poisson's ratio 

frequency 

slope of elastic axis 


Subscripts: 


ij 

Ji 


quantity in ith rib at Jth spar 
quantity in jth spar at ith rib 
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in quantity in ith rib at center of nth bay * 

quantity in jth spar at center. of mth bay 
ran quantity on spanwise cross section at center of cell mn 

nm quantity on chordwise cross section at center of cell mn 

i value in ith rib 

J value in jth spar 

K value in rib along line of support 

m value in mth bay along spar 

n value in nth bay along rib 

x value in web of spar 

y value in web of rib 

xy quantity having same value in rib or spar 

max maximum value 


PHYSICAL ASSUMPTIONS OF THEORY 


All stresses are assumed to be in the elastic range. The thickness 
of all webs is assumed to be sufficiently small that the variation of 
stress over the thickness may be neglected. It is assumed throughout 
the analysis that chordwise cross sections of the wing have a horizontal 
axis of symmetry. Some violation of this limitation can be tolerated in 
most cases without reducing significantly the accuracy of the solutions. 

Two simplifying approximations must be introduced as follows: 

(1) Since the depth of a wing varies in all directions, the shear 
flows in the skin on any cross section will have vertical components. 

The contributions of these vertical components to the shears on the 
cross sections are assumed to be negligible.. 

(2) Because of the variable depth of the wing, the normal stresses 
In the skin will have vertical components. The contributions which these 
vertical components make to the shears on the cross sections are assumed 
to be negligible. 
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As a result of the above assumptions the difference equations for 
the structure may he derived by considering a wing of constant depth 
without any loss of generality. The wing will also be assumed to have 
a constant chord as shown in figure l(a) . The wing may be supported 
at the face of the fuselage as shown. For purposes of analysis the 
original, structure may be replaced by an idealized structure as indi- 
cated in figure l(b) . All of the normal stresses are assumed to be 
carried in flanges located at the junction. of the ribs or spars and the 
skin. All connecting webs carry only shear flows. 


' EFFECT OF AXIS OF SYMMETRY OF CROSS SECTIONS 


When a wing is acted upon by vertical loads or torsional loads and 
the cross sections have a horizontal axis of symmetry, it can be shown 
that the normal stresses have an antisymmetric distribution about the 
axis of symmetry. Hence the normal stresses may be considered to form 
bending moments just as in an ordinary plate. Similarly the shear flows 
in the skin may be shown to be antisymmetrical and, thus, may be consid- 
ered to form twisting moments corresponding to ordinary plate theory. 
Since the internal stresses due to vertical or torsional loading may be 
expressed in terms of vertical shears, bending moments, and twisting 
moments, the equations which govern the structure take the form of plate 
equations. 


VERTICAL EQUILIBRIUM OF AN ELEMENT 


The first step in the development of the theory is to determine the 
relation between the internal shears and the applied load. The concen- 
trated load acts at the intersection of the ith rib and the jth spar. 

A small portion of the wing surrounding this point may be removed for 
analysis. If the skin and idealized flanges are removed from the ele- 
ment, the remaining webs will appear as shown in figure 2. Only the 
vertical forces which act on the element are shown. In the idealized 
structure the shear in a rib or spar will be constant over the width 
of a cell, or bay. The first subscript indicates the member in which 
the shear acts and the second subscript indicates the bay in which the 
shear occurs. From figure 2 the equation of equilibrium is seen to be 
as follows: 


V j(m+1) " v jm + v i(n+l) " v in + p ij = 0 t 1 ) 


Differencing symbols will now be introduced. The symbol A^ indi- 
cates the jump in a function across the ith rib. The symbol Aj indicates 
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the jump in a function across the jth spar. An application of these 
operators to the shears yields the following formulas: 

Ai^jm = Vj(m+l) - Vj m (2a) 

A j v in = Vi( n +i) - V in (2b) 

Substituting equations (2) into equation (l) gives 

AiV Jm + AjV ln = -Pij (3) 

This equation is a first-order difference equation. 


EQUILIBRIUM OF MOMENTS ON AN ELEMENT 


The next relation to be considered is that which holds between 
moments and shears. For this purpose an element of the wing must be 
considered as shown in figure 3 • This element extends along the Jth 
spar over the length Xq of the mth bay. The only forces shown on 
the drawing are those which form couples having axes that are perpen- 
dicular to the spar. From figure 3 it is seen that the equation of 
equilibrium for moments is as follows: 


Mj(i+l) - + Tjujj - T m ( n _]_) - XmVjjn - 0 (4) 

The twisting moment T mn is the total twisting moment over the width 
of a cell rather than the twisting moment per unit of width which occurs 
in ordinary plate theory. The bending moment Mji is the moment in the 

jth spar of the idealized structure at the point of intersection with 
the ith rib. 

One may consider an element corresponding to the illustration of 
figure 3 which is oriented in the perpendicular direction. This element 
would contain the nth bay of a rib of length X n with the adjoining 
skin. Equilibrium of moments for couples having axes parallel to the 
x-axis gives the following equation: 

^i(j+l) ” ^ij - Tqjjj + m— l) “ ^n^in = 0 (5) 

Additional differencing symbols may now be introduced. The symbol 
Ajn indicates the increase' in a function across the mth bay of a spar. 
The symbol A n indicates the increase in a function across the nth bay 
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of a rib. These operators may be applied to the bending moments to 
obtain the following formulas: 

AmMji = Mj(i+i) - Mji (6a) 

A n Mij = Mi( j+i) - Mij (6b) 

The operators Ai and Aj, which have been previously defined, may be 
applied to the twisting moments to obtain 

A-jTnni = T rim “ ^n(m-l) (7 a ) 

^j^mn = ^mn “ ^m(n-l) 

Equations (6) and (7) may be substituted into equations (4) and (5) to 
obtain the following first-order difference equations: 

AinMji + AjTnm = XmVjm (8a) 

AnMij - AiTnm = X n Vi n (8b) 

Since the shear flow in the skin on two perpendicular planes must 
be equal, there exists a simple relation between T mn and T nm which 
is expressed as follows: 


^mn _ _ ^nm 


(9) 


RELATION BETWEEN BENDING MOMENTS AND ROTATIONS OF NORMALS 


The next step in the development of the theory is to form relations 
between bending moments and rotations of the normals to the middle sur- 
face. .This requires a consideration of the deformation of a rib or spar 
which takes place when shearing strains are significant. In figure 4(a) 
there Is shown a segment of a spar acted upon by shears and bending 
moments. A solid vertical line has been drawn normal to the elastic 
axis before distortion. 

The spar is shown in its distorted position In figure 4(b) . The 
original normal has rotated through an angle |3 X and makes this angle 
with the vertical dashed line cd. The dashed line ab has been drawn 
perpendicular to the elastic axis in its deflected position. The original 
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normal makes an angle 7 X with ab where y x is the shearing strain 
in the web of the spar. The angle between ab and cd is equal to the 
slope of the elastic curve. 


The normal strains in the idealized flanges are proportional- to 
the rate of change of rotation of the original normals in the spar. 
Because of the Poisson ratio effect in the skin the rotation of the 
normals in the ribs will also contribute to the normal strains in the 
spars. Similar considerations can be given to a rib. This permits the 
writing of bending moments in terms of rotations of the normals in the 
ribs and spars. 


The following equations may be written for the bending moments in 
the ith rib and the jth spar, respectively, at their intersection point: 


M ij 


= -D 


ij 


/dp y \ /8p x \ 

Ur*) -DiJ’M) 

W/i i 


ij 


(10a) 


M 


Mjl = -MaTAj 


J^y\ 




(10b) 




The coefficient Dij is the bending stiffness of the rib including 

the vertical web and the skin extending to the center of each adjacent 
cell. If I is the moment of inertia of the skin per unit of width 
and I r is the moment of inertia of the vertical web of the rib, the 
formula for the bending stiffness of the rib becomes 


D ij = H — -p + Elr 
1 - H 


(11a) 


If I 6 is the moment of inertia of the vertical web of the spar, the 
total bending stiffness of- the spar becomes 


D Ji - + Ms 

l - n 


(lib) 
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The Poisson ratio effect occurs only in the skin and hence the coeffi- 
cients Dij’ and Dji' are defined hy the following equations: 




T 



(12a) 


Dji‘ 



(12b) 


For a nonuniform wing the values of I, I r , and I s may he computed 

at the point of intersection of the rib and spar. For better accuracy 
I may be computed as an average value over an appropriate area and 
I r and Is may be computed as average values over an appropriate length. 


The derivatives which are contained in equations (10) must be 
replaced by finite differences. If is the rotation of a normal 

in the ith rib at the center of the nth bay, the approximation to the 
first derivative is given by 


&L = ^ [pin ■ pi(n - 13 

The corresponding equation for the spar is given by 


(13a) 



(13b) 


If the appropriate differencing symbols are introduced, equations ( 13) 
become 


(ilrt Ajl3ln (14a) 

(sr)ij ' M (lto) 
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The finite differences which are given in equations (l^) mpy 
he substituted into equations (10) to obtain the following equations 
relating bending moments and the rotations of the normals: 


u 11 

T^AjPin + 


Dlj' 

X± 


AiPjm. - -Mij 


(15a) 



(15b) 


RELATION BETWEEN TWISTING MOMENTS AND ROTATIONS OF NORMALS 


A relationship must now be derived between the twisting moments 
and the rotations of the normals to the middle surface. The twisting 
moments may be related to average values of the shearing strains in 
the skin panels. Simple relations may be found between these average 
strains and the rotations of the normals in the ribs and spars. These 
relations can be developed readily by reference to the idealized structure. 

In figure 5 there is shown a panel of skin enclosed between two 
adjacent ribs and two adjacent spars. The positxon which the panel 
takes when the structure is loaded is shown in dashed lines. It will 
be assumed that the rotations of the normals in the ribs and spars are 
defined at midpoints which are indicated as A, B, C, and D in figure 5* 

The new positions which these four points take in the loaded structure 
are also shown. The components of displacement which these points 
take, in the direction of the members on which they are located, are 
indicated as a, b, c, and d. 

The average shearing strain in the panel may be assumed to be 
equal to the change in angle made by the two center lines. These center 
lines rotate through the angles and a n . Hence the shearing strain 

in the panel is given by 

7mn = <% + °n (l6) 

These angles can be expressed in terms of the displacements and panel _ 
dimensions as follows: 


c - a 
lm 


ctjn = 


(17a) 



NACA TN 2786 


11 


«n = (17b) 

A -n 


The displacements can he expressed in terms of the rotations of normals 
as follows: 


a - ” 2 (l8a) 

b = " 5 - £(j+l)m ( l 8 h) 

c = - §• 3(i+l)n (18c) 

d = - | Pjm (l8d) 

Substituting equations (l8) into equations (17) gives 

^ = rj^[P(i+l)n “ Pin} ( 19a) 

an = ^^I^(«3 +1 ) m ” (19b) 

Introducing the appropriate differencing symbols gives 

^ AmPin (20a) 

eA. m 

ttn = 2^ AnPjm ■' (20b) 


Equations (20) may be substituted into equation (l 6) to obtain the 
following formula for the- shearing strain: 




( 21 ) 
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The twisting moments must now be expressed in terms of the shearing 
strains. The shearing stress is obtained from the strain by multiplying 
by the shearing modulus. The shear flow is then obtained by multiplying 
by the skin thickness. 


Qmn = ' ^ 

(22) 

The twisting moment per unit of width is obtained by multiplying by the 
depth of the wing. The total twisting moments for the cell are then 
obtained by multiplying by the length or width of the cell. The formulas 
thus obtained are 

AaPin + A n 3 - -Tmn 

(23a) 

2&S. Z^aSin + A n 3 jm = T n m 

- 

(23b) 

The stiffness coefficients which appear in equations (23) may be 
from the following formulas: 

computed 

Dmn = ^m = ^m^I 

(24a) 

Dnm = ^n ^ 

(24b) 


RELATION BETWEEN DEFLECTIONS AND ROTATIONS OF NORMALS 


It now remains to express the deflections in terms of the rotations 
of the normals by finite difference equations. Returning to figure 4 
and considering similar conditions in the chordwise direction, the fol- 
lowing two equations can be written: 


~ = Px + 7 X (25a) 


8w 

8y 


+ h 


(25b) 
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The shearing strains are average values over the vertical webs of the 
ribs or spars. In the idealized structure these strains are constant 
over the shear webs. 


The shearing strains are related to the shears by the following 
formulas : 



(26a) 


„ZiL..Zz 

GA-r Cv 


(26b) 


In equations (26) Ag and Ar are the areas of vertical webs of the 
spars and ribs, respectively. Equations (26) may be substituted into 
equations (25) to obtain 


— = + 5 s - 

Bx 


(27a) 


8w _ a ,^ 2 L 

w ■ % t 


( 27 b) 


It is now necessary to replace equations ( 27 ) by difference equa- 
tions. The quantities on the right-hand side may be regarded as values 
at the center of a bay or they may be. regarded as average values over 
a bay. The left-hand side may be replaced by finite differences to 
obtain the following equations: 





(28a) 

Cjm 

V in 

^in 

(28b) 

be computed at 

the center 


of the bay although better accuracy will be obtained by computing average 
values over the bays. Introducing appropriate differencing symbols into 
equations (28) gives - 

■fe = P ^ ra + 


(29a) 
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r* An w 1(j = Pin + (29b) 

A n u in 

The complete solution is now defined by nine first-order difference 
equations. 'The nine unknowns consist of six internal force quantities 
(shears, bending moments, and twisting moments) and three displacement 
quantities (deflections and rotations of normals). The first-order 
equations may be combined to form equations of higher order. This would 
be necessary if the solution were to be computed on a digital computer. 
However, when the solution is to be computed on an analog machine the 
equations may be left in their present form. Coupled circuits may be 
designed in which the voltages and currents are related by first-order 
difference equations which have the same mathematical form as the struc- 
tural equations. From these circuits one may read any one of the nine 
unknown quantities independently. 


VIBRATION MODES 


The natural vibration modes and frequencies of the wing may be 
easily determined by replacing the distributed mass by an equivalent 
set of concentrated masses. These masses must be considered as being 
located at the intersection points of the ribs and spars. Using 
D'Alembert's principle and assuming sinusoidal motion at a frequency 
cn, the concentrated load becomes 

p iJ " ^“i^ij (3°) 

Equation (30) may be substituted into equation (3) to obtain 

A i v ,)m + A j v in - (31) 

Equation (31) must be used in place of equation (3) for vibration anal- 
ysis. All other equations remain unchanged. 

BOUNDARY CONDITIONS 


When the computations are to be made on an analog computer the 
boundary conditions on free edges may be stated directly in terms of 
force quantities. In stating the boundary conditions for quantities 
which are defined at the center of a bay or a cell, it is convenient 
to introduce a row of imaginary cells surrounding the plan form of the 
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wing as shown in figure 6 . For convenience it will be assumed that the 
plan form has double symmetry and the boundary conditions will be stated 
for the first quadrant. 

If n max is the maximum value of n on the plan form, the exterior 
row of imaginary cells in front of the leading edge may be numbered 
n max + 1 as shown in the figure. Similarly the row of cells beyond 
the wing tip is numbered + 1. The rib at the tip is numbered 

i mHY (maximum value of i) and the spar at the leading edge is num- 
bered . j max (maximum value of j) . The rib at the line of support is 
numbered K. In terms of the above-defined symbols the boundary condi- 
tions can be expressed as follows: 

Line of support: 

= 0 

Wing tip: 

V «3(“max + 1) = 0 
M «3( 3-max) ~ ® 

Tn(m mfl.x + l) = 0 

Leading edge: 

^(Dthhv + 1) = 0 

Mi (jmax) ° 

Tm (%iax + 1) = 0 

DESIGN OF ANALOGOUS CIRCUIT 


(32) 

(33a) 

(33b) 

(33c) 

(34a) 

(34b) 

(34c) 


An explanation of the design of the analogous circuit will be given 
in a manner corresponding to the treatment of beam analogies which is 
contained in reference 1. Since the design procedure has been described 
in complete detail in reference 1, the present design will be stated 
more briefly. Since the structural equations have the same form as 
plate equations the analogous circuits will be approximately the same 
as those given by MacNeal (reference 2) . The only difference will be 
due to the allowance for the effect of shearing strains upon deflections. 
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In reference 1 the structural quantities and the electrical quan- 
tities were represented- by completely separate sets of symbols. In the 
present paper the structural symbols will also be used as electrical 
symbols. Thus it will not be necessary to write additional equations 
for the electrical circuit. For convenience in the discussion of the 
circuit the structural equations which govern the problem are repeated 
in a group as follows: 


A i V jm + A j v in = - p ij 

(3) 

AflMji + 

(8a) 

A n M ij " ^i T nm = \F±n 

(8b) 

^mn _ ^nm 

(9) 

D ii D il‘ 

if Vin + Tr A l e > ■ -«ij 

(15a) 

ff A jPin + A iPjm = -Mji 

(15b) 

2m AqPjh, + 2aa ^ ln = . Tmn 

A. n 

(23a) 

if + 

(23b) 

5*^ = + 

( 29 a) 

f Ai»iJ - Pin + Jt 1 
*n Cin 

( 29 b) 


It is seen that equation (3) expresses a relation between force 
quantities and hence is represented electrically by an application of 
Kirchhoff's law to a nodal point. This nodal point is shown in fig- 
ure 7 where the structural equation is also given. The currents are 



NACA TN 2786 


17 


designated with structural symbols so that the structural equation may 
be interpreted as an electrical equation. 

The next equation to be considered is equation (8a) . This again 
calls for an application of Kirchhoff 's law at a node and is illus- 
trated in figure 8. In this figure it will be noted that the "loading 
current" proceeds from a transformer through a resistance. 

This resistance is required to allow for the effect of shearing strains 
upon the deflections. Equation (8b) is of a similar nature and the 
analogous node is shown in figure 9* 

Three different nodal points have been illustrated. Each of these 
nodes must occur in a different planar circuit. Hence three p lan ar 
circuits must be designed and interconnected by means of transformers. 

A panel of the shear circuit is shown in figure 10. In this circuit 
the voltages represent deflections and the currents represent shears. 

The design of the other two planar circuits requires the intro- 
duction of additional circuit elements. The next structural equations 
to be considered are equations (15). For more convenient electrical 
representation these equations are solved for the differences in rota- 
tions in terms of bending moments. The result of such a solution may 
be written as follows: 

A iPjm = - ®x^ji + ^xy^ij (35a) 


A j%n = RxyMji - RyMij (35b) 

The circuit elements corresponding to equations (35®) and (35b) are 
shown in figures ll(a) and 11(b) , respectively. In order to see more 
clearly that the voltages and currents shown in figure 11 are related 
by equations (35) one may rewrite equations (35) in the following form: 


A iPjm ( R x “ E x^ M ji 

- - M U) 

(36a) 

Aj3in = - R xy|Mij - 

- (Ry - 

(36b) 


Equation (38a) may be immediately checked by an inspection of 
figure ll(a). In order to see that equation (38b) is valid one must 
use the transformer laws. These laws have been illustrated in refer- 
ence 1. One-half of a transformer is shown in figure 11(a) and the 
other half of the same transformer is shown in figure ll(b) . This 
transformer has a 1:1 turns ratio. The plus and minus signs indicate 
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the polarity (or phase) of the terminal voltages. The manner in 
which the circuit elements of figure 11 are coupled is illustrated 
in figure 12. 

The next equations to he considered are equations (23) • For 
convenience in designing the circuit these equations may he written 
in the following modified form: 

A mPin = - Tmn (37a) 

^n^jm + ^m^in = (37h) 


Equations ' (37a) and (37b) are not independent since the right-hand 
sides are related as shown hy equation (9)* A circuit element is shown 
in figure 13 (a) in which the voltages and currents are related hy equa- 
tion (37a). The associated circuit element is shown in figure 13(h). 

The manner in which the elements of figure 13 are coupled is shown in 
figure 14. A verification of equation (37a) as an electrical equation 
requires an application of the transformer law for voltages. Equa- 
tion ( 9 ) is a direct statement of the transformer law for currents. 

It is seen from these equations that the ratio X n /l m must he the 

transformer turns ratio. 

The next equations to he considered are equations (29) • These 
equations may he interpreted as electrical equations hy referring to 
figures 15 and 16 . Figure 15 shows the transformer connection between 
the shear circuit wherein voltages are deflections and the spanwise- 
bending-moment circuit wherein the voltages are rotations of the normals 
In the spars. The point at which the voltage corresponds to the slope 
of the elastic curve of the spar is indicated as Figure l 6 shows 

a similar connection between the shear circuit and the chordwise-bending- 
moment circuit. 

A panel of the circuit for spanwise bending moments is shown in 
figure 17« A similar circuit for chordwise bending moments is shown 
in figure l 8 . Complete assembled circuits are shown in reference 2. 


CONCLUDING REMARKS 


A theoretical method of analysis has been developed whereby the 
stresses and deflections may be determined for a thin supersonic wing 
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of low aspect ratio. It is assumed that chordwise cross sections of 
the wing have a horizontal axis of symmetry. This permits the reduction 
of the solution to the form of plate theory. In order to carry out the 
analysis on an analog computer the theory is expressed in terms of first- 
order difference equations. The analogous electrical circuits are also 
derived. 


California Institute of Technology 

Pasadena, Calif., February 5, 1952 
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(a) Segment before distortion. 



(b) Segment after distortion. 


Figure 4.- Distortions of a segment of a spar. 




NACA TN 2786 


25 


I 



Figure 6. - Symbols required for statement of boundary conditions 
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Figure 7-- Nodal point in shear circuit. AjV + AjV in = -P.y. 



Figure 8.- Nodal point in spanwise-bending-moment circuit. 

AjjjMji + AjT n m = 
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-*-A/V 

Pjfr-0 
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Wji- Mjj 


*jm 


(a) Spanwise bending moments, 
A i^Jm = “^x^ji + ^xy^i j ■ 



(b) Chordwise bending moments. 
A j^in = ^yMji " R y M ij* 


Figure 11.- Circuits for bending moments. 



Figure 12.- Coupling of circuits for bending moments. 
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(a) Twisting moment on spanwise ("b) Twisting moment on chordwise 

cross section. , cross section. 


Figure 13- - Circuits for twisting moments. = - — — 
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Figure l4. - Coupling of circuits for twisting moments. 
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Figure l6. - Relation between deflections and chordwise slopes. 
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Figure 18. - Panel of circuit for chordwise tending moments. 
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